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Abstract 

The observation of air showers from elementary particles with energies exceeding 10 20 eV 
poses a puzzle to the physics and astrophysics of cosmic rays which is still unresolved. 
Explaining the origin and nature of these particles is a challenge. In order to constrain 
production mechanisms and sites, one has to account for the processing of particle spectra 
by interactions with radiation backgrounds and magnetic fields on the way to the observer. 
In this paper, I report on an extensive study on the propagation of extragalactic nucleons, 
7-rays, and electrons in the energy range between 10 8 eV and 10 23 eV. I have devised an 
efficient numerical method to solve the transport equations for cosmic ray spectral evolution. 
The universal radiation background spectrum in the energy range between ~ 10 _9 eV and 
~ 1 eV is considered in the numerical code, including the diffuse radio background, the 
cosmic microwave background, and the infrared/optical background, as well as a possible 
extragalactic magnetic field. I apply the code to compute the particle spectra predicted by 
various models of ultrahigh energy cosmic ray origin. A comparison with the observed fluxes, 
especially the diffuse 7-ray background in several energy ranges, allows one to constrain 
certain classes of models. I conclude that scenarios which attribute the highest energy 
cosmic rays to Grand Unification Scale physics or to cosmological Gamma Ray Bursts are 
viable at the present time. 

PACS numbers: 98.70. Sa, 98.70.Rz, 98.80. Cq 



1 Introduction 



Shortly after the cosmic microwave background (CMB) was discovered H it became clear 
that this universal radiation field has profound implications for the astrophysics of ultrahigh 
energy cosmic rays (UHE CR) of energies above 10 18 eV. For nucleons the most profound 
effect is photoproduction of pions on the CMB. Known as the Greisen-Zatsepin-Kuz'min 
(GZK) "cutoff" 0, [J, this effect leads to a steep drop in their energy attenuation length by 
about a factor 100 at around 6 x 10 19 eV which corresponds to the threshold for this process. 
The nucleon attenuation length above this threshold is about 10 Mpc. Heavy nuclei with 
energies above about 10 19 eV are photodisintegrated in the field of the CMB within a few 
Mpc M. One of the major unresolved questions in cosmic ray physics is the existence or 
non-existence of a cutoff in the UHE CR spectrum at a few 10 19 eV which, in the case of 
extragalactic sources, could be attributed to these effects. 

Therefore, there has been renewed interest in UHE CR research since events with energies 
exceeding 10 20 eV have been detected. The Haverah Park experiment [0] reported several 
events with energies near or slightly above 10 20 eV. The Fly's Eye experiment |8], [| detected 
the world's highest energy CR event to date, with an energy ~ 3 x 10 20 eV. Near the 
arrival direction of this event the Yakutsk experiment [|1(J recorded another event of energy 
~ 1.1 x 10 20 eV. More recently, the AGASA experiment [11, 12] has also reported an event 
with energy 1.7 — 2.6 x 10 20 eV. It is currently unclear whether these events indicate a 
spectrum continuing beyond 10 20 eV without any cutoff or the existence of a cutoff followed 
by a recovery in the form of a 



: 'gap" in the spectrum 
There has been much speculation about the nature and origin of these highest energy 
cosmic rays (HECRs) ||, |l|, U, [l7|, [0 



Concerning the production mechanism one can 
distinguish between two broad classes of models: Within acceleration models, charged pri- 
maries, namely protons and heavy nuclei are accelerated to very high energies |19 , |20| in a 
"bottom-up" manner. Preferred sites are large-scale astrophysical shocks which occur for 
instance in radio galaxies ||21|| . Even there it seems barely possible to accelerate CRs to the 
required energies JI3], 122], 1231 . Recently it has also been suggested that acceleration of UHE 



CRs could be associated with cosmological gamma-ray bursts (GRBs) p4|, |26|. In the 
second class of so called "top-down" models, charged and neutral primaries are produced at 
UHEs in the first place, typically by quantum mechanical decay of supermassive elementary 
X particles related to grand unified theories (GUTs). Sources of such particles at present 
could be topological defects (TDs) left over from early universe phase transitions caused by 
the spontaneous breaking of symmetries underlying these GUTs |27l f 



29, mm 



The injection spectra in top-down models tend to be considerably harder (flatter) than in 
acceleration models. 

The particle identity of the UHE CRs is not known either. The Fly's Eye analysis || 
suggested a transition from a spectrum dominated by heavy nuclei to a predominantly light 
composition, i.e. nucleons or even 7-rays, above a few times 10 19 eV. However, this has 
not been confirmed by the AGASA experiment [Q. Although there have been claims that 
the shower profile of the highest energy Fly's Eye event may be inconsistent with a primary 
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photon |37]] or even with a proton primary f55fl , the situation is not settled because of many 
uncertainties which can affect the shower development in the atmosphere. 

Other options discussed for the nature of the HECRs include heavy nuclei and even 
neutrinos [[37J. Heavy nuclei have their own merits because they can be deflected considerably 
by the galactic magnetic field which relaxes the source direction requirements [§, O. In 
addition, for shock acceleration, heavy nuclei can be accelerated to higher terminal energies 
because of their higher charge. However, one should note that the range for heavy nuclei 
is limited to a few Mpc as mentioned above. Neutrinos, on the other hand, do not lose 
much energy over cosmological distances |39|, IQfl , but by the same token the probability for 
interacting in the atmosphere is small. Attributing the HECRs to neutrinos would therefore 
require a neutrino flux at UHEs which is much higher than the observed CR flux at the 
same energies. This poses severe constraints on the possible sources for these neutrinos pT|. 
In addition, neutrinos would be expected to give rise to predominantly deeply penetrating 
showers in the atmosphere. 

The production spectrum of UHE CRs is modified during their propagation. There are 
many studies on nucleon propagation in the literature using analytical f|3, |5| as well 

as numerical approaches [|17| , [46| , [47| , [48[| , and the propagation of heavy nuclei has also been 



considered ^L7j. This was mainly motivated by the conventional acceleration models which 
usually predict UHE CR fluxes to be dominated by these particles. However, secondary 
7-rays and neutrinos can also be produced, for example as decay products of pions created 
by interactions with various radiation backgrounds at the source or during propagation E7| . 
Under certain circumstances their flux can become comparable with the primary flux [ff9| . 
Furthermore, within TD models 7-rays are expected to dominate to begin with [Q. A study 
on 7-ray propagation in this context has been performed recently |jT] using a quantitative 
treatment on the cascade initiated by UHE photons. In my opinion, however, it suffers from 
several unrealistic assumptions with respect to the injection scenarios considered. I improve 
on their treatment of the propagation of 7-rays. Apart from that I find three reasons to 
explore UHE 7-ray propagation in more detail in this paper: 

First, due to the absence of threshold effects similar to photopion production which causes 
the GZK "cutoff" for nucleons, the 7-ray spectrum is not expected to have a break around 
10 20 eV. Furthermore, 7-rays can generate electromagnetic (EM) cascades while propagating 
rather than being absorbed right away. UHE electrons produced by pair production upscatter 
background photons and transfer most of the energy back to photons. This effect consider- 
ably increases the effective energy attenuation length of the "cascade" photons |52], |53j . At 



a few times 10 20 eV this attenuation length may be even greater than that for protons which 
drops precipitously at the threshold for photopion production. Extragalactic 7-rays could 
therefore have some potential to produce a recovery beyond the GZK "cutoff" . 

Second, in contrast to the case of nucleons, the propagation of 7-rays is presently fraught 
by certain ambiguities which are mainly due to uncertainties in the intensity of the uni- 
versal radio background and the strength and spectrum of the extragalactic magnetic field 
(EGMF). We hope that an application of the general framework presented here under dif- 
ferent assumptions for such parameters could in turn provide some insights into their actual 
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values once the UHE 7-ray flux is known to some accuracy. This would be in some anal- 
ogy to the method of using TeV 7-ray observations to constrain or detect the universal 



infrared/optical background [5^|. In previous work [52] it is shown that, depending on its 
strength, the large-scale EGMF could produce a feature in the 7-ray spectrum which might 
be observable in the future. 

Finally, the study of high energy cosmic and 7-ray propagation can place stringent con- 
straints on the nature and origin of UHE CRs. Such constraints can be obtained by com- 
puting the propagation modified spectra especially of lower energy 7-rays expected within 
a certain scenario and comparing the predictions with the observed fluxes [|49| , [5"5| , [56|. At 



UHEs there are some experimental prospects to distinguish 7-rays from other primaries in 



the future, possibly even on an event by event basis |57 [. This would allow comparing not 
only the total fluxes of UHE nucleons, heavy nuclei, and 7-rays, but also their composition 
with model predictions. 

This motivated the present comprehensive study of propagation of nucleons and 7-rays 
and its application to models which attribute UHE CRs to top-down mechanisms within 
GUT-scale physics or associate them with cosmological GRBs. I explore the energy range of 
10 8 < E < 10 23 eV. The low end is chosen such that we can draw constraints by comparing 
the propagated spectra with existing measurements of the diffuse 7-ray background around 



100 MeV |5^, |59], p0| . The high end is chosen beyond the highest CR energies ever observed 
enabling us to study top-down models. I include not only the CMB but also the diffuse radio 
background which plays a big role at the highest energies and the infrared/optical (IR/O) 
background which influences the flux at somewhat lower energies. I also include the EGMF 
as a free parameter. The propagation of nucleons is also studied with special emphasis on 
the production of secondary 7-rays, electrons, and neutrinos. 

The rest of the paper is organized as follows: In Section 2, 1 present the general ingredients 
of calculating the propagation of extragalactic 7-rays and nucleons. I discuss the role and 
nature of the low energy photon background and the EGMF, and explain in detail the 
implicit method used in solving the transport equations numerically. Section 3 is devoted 
to the treatment of the relevant interactions of 7-rays and nucleons. I compare our analysis 
with other work in Section 4. Section 5 discusses the generic forms of the injection spectra 
and the source distribution for typical top-down models and the GRB scenario. Results 
and constraints from the spectra predicted at Earth are presented in detail. In Section 6, I 
summarize the findings and discuss future prospects. 

2 Formalism 

2.1 Radiation backgrounds 

UHE CRs undergo reactions with the universal diffuse radiation backgrounds permeating 



the universe [p2 |. The most relevant among them are the CMB, and the radio and IR/O 
backgrounds. 

Photon primaries of energy E can be absorbed by pair production with a background 
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photon of energy e if (c = % — 1 throughout) 

E > E th = ^ = 2.611 x 10 11 (jL^j 1 eV . (1) 

Therefore, for a typical CMB photon (e ~ 10 -3 eV) the threshold energy Eft is ~ 3 x 10 eV, 
whereas for a typical radio photon (e < 1CT 8 eV) the threshold is > 3 x 10 19 eV, thus 
affecting UHE 7-rays. Furthermore, since the pair production cross section peaks near the 
threshold, pair production on the radio background dominates over pair production on the 
CMB in that energy range although the number density of radio background photons is much 
smaller than that of CMB photons. On the other hand, the IR/O background affects the 
lower energy photons for the same reason. The threshold for pair production on the IR/O 
background lies at about 10 12 — 10 13 eV. Similarly, the contribution of the IR/O background 
to the total photopion production rate by protons is not negligible in the lower energy range 
(£<10 18 eV). 

All these backgrounds evolve with time (i.e. redshift) by cooling with the expansion of 
the universe. However, on top of it, the radio background and the IR/O background evolve 
due to the evolution of the respective sources. The evolution of the radio background is 
tied to the evolution of the radio sources such as radio galaxies, and the evolution of the 
IR/O background to that of normal galaxies. Treating the evolution of these backgrounds 
carefully is important if we are to go back to the redshift where there existed not many of 
these sources (z > 5 — 6). The flux of an isotropic radiation background component produced 
by an ensemble of sources is given by the following relation: 

where j(e, z) is the radiation flux (in units of number per area per time per solid angle per 
energy) at redshift z, Z{ is the initial redshift when the sources begin to appear, H is the 
Hubble constant, and $(e, z) is the production spectrum of the relevant background (in units 
of number of background photons per volume per time per solid angle per energy) at redshift 
z. Throughout this paper we assume Hq = 75 km sec -1 Mpc -1 and a critical density universe 
(i.e. Qq = 1) for simplicity, but I keep Hq in the formulae to show the dependence. If we 
assume that 

^(e,z) = %(e)(l + z) 3 N c (z), (3) 

where $o( e ) is the typical intensity spectrum of an individual source and N c (z) is the co- 
moving density of the sources as a function of the redshift, Eq. (fj) may be rewritten as 

j(e, z) = ^(1 + zf Jjl + z')- 3/2 $o (e i±|) N c {z')dz> . (4) 

The background photon number density is then obtained from the relation n(e, z) = 4tt/c ■ 
j(e,z). Note that it is important to self-consistently derive the background in this way 
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rather than following the often-used approach of assuming a current background photon dis- 
tribution, n (e), and then extrapolating it back to higher redshifts via the formula n(e, z) = 
(1 + z) 3 no[(l + z)e\. While easy to implement, this approach is only valid for a truly pri- 
mordial background formed at extremely high redshifts (e.g. the CMB) and can lead to 
misleading results if one is not careful. 

The role of the IR/O background in determining the level of the cascade radiation back- 
ground below ~ 1 TeV which is the diffuse 7-ray background due to cascades by CRs was 



examined in more detail in Ref. [63]. Even for rather extreme assumptions for the IR/O 
background, the cascade background level typically does not vary by more than a factor of a 
few. Accordingly, for this paper I have chosen to adopt a simple, "middle of the road" model 
for the formation of the IR/O background (a discussion of the various possibilities can be 
found in Refs. |63|, |65|). I assumed that the dominant contribution to the IR/O background 
comes from ordinary galaxies which formed early in the universe, at Zi ~ 5. The typical 
galaxy was assumed to have a spectrum like that of the 5 Gyr disk galaxy spectrum shown 



in Fig. 4 of Ref. |^6| , which has a component peaking at ~ 1 fim in wavelength due to direct 
emission from stars and a second component peaking at ~ 100 /im due to reprocessing of the 
starlight by interstellar dust. The combined number and luminosity evolution of the galaxies 
was taken to go as (1 + z) 7 , i.e. most of the background was produced in a strong, initial 
burst of star formation in the galaxies, and the intensity of their emission was adjusted to 
give an optical background density today of n opt ~ 2 x 10~ 3 cm -3 . 

For the present diffuse extragalactic radio background spectrum I use the estimate given 
in Ref. |67] (see also Refs. |62|, ^8|). This spectrum can be parametrized by a power law with 
a lower frequency cutoff for which I use f c = 2 MHz. One can also estimate the contribution 
to this background in the power law regime caused by radio galaxies. I did that by inserting 
the injection flux $(e, z) resulting from the radio luminosity function given by Eq. (7) of 



Ref. |69[ into Eq. (|2]). The intensities resulting at z = are within a factor ~ 2 of the 
estimate given in Ref. [p7 ] . I adopt the functional redshift dependence for the power law 
regime following from this calculation and normalize it to the present intensities given in 
Ref. |67| . In addition, I assume a redshift-independent lower frequency cutoff at f c = 2 MHz. 

The combined radiation spectrum at z = used in this paper is presented in Fig. 1. In 
Fig. 2 I plot N c (z), i.e. the effective comoving densities of radio and IR/O sources whose 
luminosities are normalized at z = 0, as functions of redshift. 



2.2 Extragalactic Magnetic Field 

The long range EGMF affects the propagation of CR particles via synchrotron radiation and 
deflection (or even diffusion). 

Synchrotron radiation is much more straightforward to consider than deflection. The 
synchrotron loss rate for a charged particle with mass m, energy E, and charge qe (e is the 
electron charge) subject to a magnetic field of strength B is given by |61] 

dE * fV£M7*V. (5 ) 



dt 3 87r V m ) \m. 
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where <jt is the Thomson cross section, and m e is the electron mass. Here, the average over 
random magnetic field orientations was taken. Synchrotron loss influences the electronic 



component of the cascade most strongly in the UHE regime [52]. On the other hand, at a 
given energy the synchrotron loss rate for protons is much smaller than that for electrons 
because the loss rate is proportional to m -4 . Thus, for protons synchrotron loss is completely 
negligible for the energies, magnetic field values, and distances I consider in this paper. 
The relevant synchrotron power spectrum radiated by the electrons is given by |H] 



dP ^/3e 3 B 

where 



dE 1 2ix m t 



G(E,/E C ) , (6) 



G(x) EE X y/l - {x/£)*K B/3 {Z)dS , (7) 

J x 

and the critical energy E c is defined as 



3eB ( E P \ _ _ „ „ 14 / E P \ ( B 



E c ee ^— — ( — — ) ~ 2.2 x 10 i4 ( ) ( ) eV . 



2m P \mj ' V10 21 eV/ V10^ 9 G 



The power spectrum peaks at E 1 ~ 0.23E C . The number spectrum, which is obtained by 
dividing Eq. (||) by the photon energy E 7 , is a monotonically decreasing function of energy. 

Deflection is another important factor when dealing with the propagation problem in 
general |52j| . The straight line propagation (SLP) approximation which treats the motion 
of CR particles in one dimension fails if the effect of the deflection becomes large. The 
gyroradius of a charged particle with charge qe and momentum p (energy E) is given by 

where B± is the field component perpendicular to the particle's motion. Note that the 
EGMF deflects protons and electrons by the same amount at a given energy once they are 
relativistic. If the gyroradius of a charged particle is considerably longer than the source 
distance, the effect of the deflection is practically negligible. On the other hand, if the gyro- 
radius is comparable or shorter than the source distance, the deflection may not be neglected 
and one now has to keep track of the transversal motion, which makes the problem much 
more complicated. However, if the sources are distributed homogenously and isotropically 
throughout the universe, then the influence of the deflection on the shape of the spectrum 
becomes small. Although this is a purely mathematical model, it is a good approximation 
for many realistic situations. On the other hand, if one considers the CR flux from a single 
source, deflection becomes important. 

In the case of an EM cascade the propagating particle basically alternates between a 
photon and an electron, and only electrons are affected by the EGMF. Therefore, the effective 
gyroradius of a cascade photon W g can be expressed as 

iq - K (i + ^) (io) 
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where R e is the electron gyroradius, L 7 is the photon interaction length, and L e is the 
energy loss length for the electron. If the effective gyroradius is considerably shorter than 
the source distance, the real spectrum would be very different from what one obtains by 
using the SLP assumption, and below the energy where the gyroradius is comparable to the 
source distance the flux is expected to be heavily suppressed. This point has been ignored 
in most of the work on CR propagation j|6], [47], ^TJ. Fig. 3 illustrates the gyroradii and the 
synchrotron loss rates of electrons for various strengths of the EGMF. In this paper, the 
strength of the EGMF is assumed as a free parameter between 10~ 12 G and 10~ 9 G as in 



Ref. [52 



2.3 Transport Equations 

I adopt a transport equation scheme to solve the propagation problem. Since we have an 
EM cascade ensuing, it is often inadequate to use the simple continuous energy loss (CEL) 
approximation which neglects non-leading particles. In addition, since particle numbers grow 
fast with time, using a full-blown Monte Carlo calculation will require excessive computing 
time. In this problem, using the transport equation approach is very economical in terms of 
computing time as well as sufficiently accurate. Previous work done by Protheroe & Johnson 
|5Tfl uses a mixture of transport equations and Monte Carlo techniques. 

A sample transport equation for electrons which includes pair production (PP) and in- 
verse Compton scattering (ICS) can be written as follows: 

d . f f 1 — /3 e /l 



dt 



N e (E e ,t) = -N e (E e ,t)Jden(e)Jdfi^^a lC s(E e ,e,fi)+ (11) 
J dE' e N e (E>,t) J den(e) J d / ~ ^ ^ (E e ; E' e , e, fj) + 
J dE y N y {E^t) J den(e) J dy ~ ^ °™ (E e ; E y , e, fj) + Q{E e ,t), 

where N e (E e ,t) is the (differential) number density of electrons at energy E e at time t, n(e) 
is the number density of background photons at energy e, Q(E e , t) is an external source term 
for electrons at energy E e and time t, [i is the interaction angle between the CR electron 
and the background photon (/i = — 1 for a head-on collision), and f3 e is the velocity of the 
CR electron. The terms describe the loss of electrons due to ICS, the influx of electrons 
scattered into the energy range due to ICS, the influx of electrons produced due to PP by 
photons, and the external injection. The factor (1 — f3 e /j,)/2 is the flux factor. I define the 
angle averaged cross sections R(E, e) and P(E'; E, e) as 

R(E,e) = ( d/j, 1 ~ ^ cr(E,e,/x), (12) 



and 

P{E'- E,e) = J d^^-^^(E'; E, e, y). (13) 
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Then Eq. flTT| ) is rewritten as 



j t N e (E e ,t) = -N e (E e ,t) J den(e)R ICS (E e ,e) + J dE' e N e (E' e ,t) J den(e)P eiICS (E e ; E' e , e) 
+ J dE^E^t) J den(e)P ejPP (E e ;E 1 ,e)+Q(E e ,t). (14) 

In order to solve this differential equation numerically, we first bin the energies of the CR 
electrons, CR photons, and background photons. We divide each decade of energy into 20 
equidistant logarithmic bins and call the central value of the z-th bin Ei and boundary values 
£i_i/2 and E i+ i/ 2 . And we replace the continuous integrals by finite sums, and integrate 
Eq. (|np over one CR energy bin. Then we get 

j i j 

^E^AeM^P^ + Q", (15) 

i j 

where N* = dE N(E,t), R k * = R(E k , ej ), P^ k = /j^ 2 dE P{E; E^), Q i = 

Se^I'/I dEQ(E,t), and Ae,- = e m/2 - ti-i/2- 

I adopt a first order implicit scheme to solve this difference equation ([15]); i.e. 



i + EyAe^O^s-P^s) 



where iV* is a solution advanced by a timestep At from iV* . Eq. (0) can be understood as 
follows: the second term in the denominator corresponds to the net loss of the particles from 
bin k. The second term in the numerator corresponds to the net scattered particle influx into 
bin k due to scattering from other bins which conserves the particle species. The third and 
the last terms in the numerator describe the influx of the particles either due to production 
by different particles or due to external injection. Various different interactions can be 
included in this main equation according to the general scheme laid out above, including 
certain energy losses which can be treated as continuous such as synchrotron radiation for 
which I use a simple first order upwind scheme. 

The implicit method has the advantage that the solution converges for arbitrary size of 
the timestep we take. Therefore, we are allowed to use a bigger timestep than is allowed by 
an explicit Euler scheme. However, to ensure the desired accuracy, we need to optimize the 



stepsize for a given problem by trial and error. See Ref. [7D| for a more detailed discussion 
of this implicit method. 

It is important to monitor conservation of particle numbers and total energy in order to 
obtain reliable results. For example, for ICS the coefficients should satisfy 

R kj = = ]T P^ k (number) (17) 
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and 

(E k + efiR* = £ EiP? k + £ (energy). (18) 

i i 

It is sometimes necessary to adjust the coefficients in order to obey these relations. This 
stabilizes the calculation against growing errors due to discretization of variables. 

In Eq. (|IED the coefficient matrices are multiplied with and summed over the background 
photon spectrum vector and/or the CR spectrum vector at a given redshift in order to 
advance the solution. This procedure has the advantage that one can deal with an arbitrary 
evolution of the radiation backgrounds in time, which is important in this problem. If one 
would integrate the background spectrum into the coefficients beforehand, it would become 
extremely difficult and time-consuming to handle an arbitrary background evolution because 
one would have to recalculate the coefficients at each redshift. 

In Ref. [51 , it was assumed that all radiation backgrounds exhibit a trivial evolution by 
redshifting. This allowed them to adopt a matrix doubling method fFI| for the propagation 
calculations. However, in case of a more realistic background evolution, matrix doubling is 
almost impossible. In contrast, our approach is always guaranteed to work efficiently and is 
sufficiently accurate in the more general case. 

Eq. (|T6| ) is then solved iteratively by inserting the initial values for N l ' J s on the right 
hand side and re-inserting the new values until a convergence is achieved. Since there are 
four main particle species (nucleons, photons, electrons, and neutrinos), one should converge 
all spectra simultaneously. However, it is economical and equally valid to converge each 
particle spectrum separately while holding the others fixed, and repeat this whole procedure 
until all spectra converge. 

In addition, we account for redshifting by performing the operation N a (E,z) — > [1 + 
Az/(1 + z)]~ 2 N a (E[l + Az/(1 + z)], z) for each particle species a after a step Az in redshift. 
Here I match Az conveniently with the logarithmic energy bin size, log 10 [l + Az/(1 + z)\ = 
\og 10 (Ei/Ei-i) = 1/20, which corresponds to the transformation N l a (z) — > [1 + Az/(1 + 
z)]-*Nl+\z). 

The numerical code is a combination and extension of the codes developed by Coppi & 
Kdnigl J70J, and by Lee & Sigl @. 



3 Interactions of Relativistic Nucleons and 7- Rays 

In this section I discuss various relevant interactions and their cross sections from which the 
coefficients R and P used in the transport equations are calculated. 

3.1 "Cascade" Photons 

Pair production (PP) and inverse Compton scattering (ICS) are the two main processes that 
drive the EM cascade. First, let us define the inelasticity which is the fraction of the energy 
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that is transferred from the scattering particle to the scattered (or produced). It is given by 



1700 = 1 — / ^Ae',s), (19) 

cr tot (s) J de' 

where s is the squared center of mass (CM) energy and e' is the energy of the recoiling 
(leading) particle in units of the initial particle energy. 

In the extreme Klein-Nishina limit where s ^> m? e) either the electron or the positron 
produced in a pair production event typically carries almost all of the initial total energy. 
The produced electron (positron) then undergoes ICS, and the inelasticity for ICS in this 
high energy limit is more than 90 %. Therefore, e~(e + ) loses most of its energy and the 
background photon is upscattered with almost all of the initial energy of the UHE photon. 
This cycle of the "cascade" photon is responsible for slowing down the energy attenuation of 
the leading particle. In some previous work it was incorrectly claimed that the UHE photons 
lose energy very fast based on the fact that the mean free path of the UHE photon is fairly 
short, but this sequence of PP and ICS makes the actual energy attenuation much slower. 
In addition, the contribution of non-leading particles to the flux which are neglected in the 
CEL approximation can be substantial for cascades which are not fully developed. This will 
be important in some of the applications considered in this paper. 

If the EGMF is present, however, the above scenario changes somewhat. In the energy 
range where the synchrotron loss rate for the electrons is greater than the ICS rate, the 
development of the EM cascade is heavily suppressed. Its penetration depth is basically 
reduced to the photon mean free path in this energy regime. A more detailed discussion is 
found in Ref. 



3.1.1 Pair Production 

The total cross section for PP (775 — > e~e + ) is well-known, and is given by 



o-pp = cr T • _ (1 - (] ) 
lb 



(3-/3 4 ) l n l±|-2/3(2-/3 2 ) 



(20) 



where /3 = (1 — Am 2 e / s) l l 2 is the velocity of the outgoing electron in the CM frame. In order 
to calculate the differential cross section, I adopt the simplifying approximation where the 
dependence of the cross section on the azimuthal angle of the outgoing particle in the CR 
frame is ignored. Since the terms that depend on the azimuthal angle is smaller by more 
than 10~ n than the leading order terms, this approximation is very accurate for practical 
purposes. The differential cross section for a photon of energy E 1 to produce an electron of 
energy E' e is then given by J72], j73| : 



da 



pp 



dEL 



EL 



Ey-E> 



+ 



EL 



+E ^-^(k + E^) 



(21) 
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where the range is restricted to (l — /3)/2 < E'JE 1 < (l+/3)/2. The differential cross section 
with respect to the positron energy is identical due to symmetry. 



3.1.2 Double Pair Production 

Double pair production (DPP; 775 — ► e~e + e~e + ) is a higher order QED process that affects 
the UHE photons. It is known that the DPP total cross section is a sharply rising function 
of s at the threshold and approaches the asymptotic value quickly at cr(oo) ~ 6.45 fib [f74"| . 
For interactions with the microwave background, the DPP rate begins to dominate over the 
PP rate above ~ 10 21 eV. If we take the contribution of the radio background into account, 
this energy goes up somewhat. 

The differential cross sections of DPP may be obtained through second order QED cal- 
culations, but it is extremely involved, and I could not find a suitable reference in which the 
differential cross section is calculated. In addition, since it is still a small sized effect, I think 
that introducing a reasonable assumption about the differential cross section is adequate for 
our purpose. Therefore, I use the assumption where one pair of the two carries all the initial 
energy and two particles in the pair share the energy equally. I presume that this assumption 
does not change the calculations in a significant way. 

In Fig. 4(a) we plot a(s) and a(s)rj(s) of which the latter is proportional to the fractional 
energy loss rate of the leading particle, for PP and DPP. 



3.1.3 Inverse Compton Scattering 

The total cross section for ICS (e7b — > c-f) is given by the well-known Klein-Nishina formula: 

2 



ClCS = C"T • 



where j3 



3m 2 e 
ml) /{a 



WW)[ 2 + V ,-f- - -1(2 - 3? - fi) lni±| 



(22) 



is the velocity of the outgoing electron in the center of mass 



frame. Most part of the energy range of interest is in the extreme Klein-Nishina regime, but 
nonetheless I use the exact formula. 

The differential cross section for an electron of energy E e to produce an electron of energy 
E' e is then given by W% [7j|: 



da 



ics _ 3mg 1 1 + 
dE' e ° T ' 8 s E e (3 



E' E. 2(1-/3) 



— + — + 

E„ EL 







p 2 



• (23) 



where the range is restricted to (1 — (3) / (1 + (3) < E' e /E e < 1. The differential cross section 
with respect to the energy E' of the outgoing photon is obtained by substituting E e — E' 
for E' e in Eq. 



3.1.4 Triplet Pair Production 



Triplet pair production (TPP; c-fb ee e + ) is a rather significant contribution to the 
interactions of UHE electrons. This process is discussed in detail in Refs. |70L WR [76], [77 1. 
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Although the total cross section for TPP on CMB photons becomes comparable to the ICS 
cross section already at ~ 10 17 eV, the actual energy attenuation is not important until 
much higher energies because the inelasticity is very small (< 10~ 3 ). Nonetheless, it is fairly 
efficient in channelling the energy content to lower energies, and may not be ignored. 

I use the formulation given by JF6| in calculating the total cross section, and the detailed 
expressions are given in Appendix A. The total cross section of TPP increases asymptotically 
logarithmically with s: 



3a 

(j TPP = o T — 

87T 



28, s 218 
— In — 

9 m 27 



s > mi), (24) 



where a is the fine structure constant. 

While it is possible to calculate the differential cross sections numerically using the ex- 
pressions given in Ref. JF5, [7^, |77| , it is extremely time-consuming because it involves multi- 



dimensional integrations of very complicated functions. Furthermore, some of the variables 
introduced there become very large or very small, and hence create problems with the fi- 
nite computing precision. The detailed behavior of the TPP cross sections near threshold is 
unimportant since TPP is dominated by ICS in this energy regime. Thus, it will suffice to 
use a simple and efficient approximation that works very well for the region away from the 
threshold. First, I make note of the fact that the differential cross section with respect to the 
energy of one of the particles of the produced pair tends to da/dE' oc E'^ 7 ^ 4 for s ^> m 2 p [75 



Furthermore, in the same regime the inelasticity for TPP can be well approximated by 173 

r)(s) ~ 1.768(s/m 2 )" 3/4 . (25) 

I then make the assumption that the differential rate P{E'; E, e) for the produced particle 
with energy E' and for the incoming electron with energy E and the incoming background 
photon with energy e is given as a power law with spectral index 5: 

P(E'; E, e) = R(E, e)C(E, e)E'- 6 , (26) 

where C(E,e) is a normalization factor. Then using the requirement that the integrated 
differential rates must be the same as the total rate and energy conservation [i.e. the ana- 
logues of Eqs. fli"7|) and fli"8|) 1, one can uniquely determine the coefficient C(E, e) and 5. The 
spectral index 5 approaches 7/4 for large s. 

For the recoiling electron, on the other hand, I may assume continuous energy loss whose 
rate is given by 

— ~-E J den{e)R{E,e)7]{s). (27) 

The importance of TPP again depends on the presence and the strength of the EGMF. 
If the EGMF is stronger than about 10 -12 G, then TPP energy loss is dominated by syn- 
chrotron cooling, and it is no longer very important. Since various arguments and indirect 
measurements of the EGMF @ suggest that EGMF is at least 1(T 12 G, TPP may not play 



a big role in the propagation of UHE photons. However, in the absence of the EGMF, the 
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contribution of TPP to the energy attenuation of electrons and photons is comparable to or 
even greater than ICS above ~ 10 22 eV and thus may not be ignored. 

In Fig. 4(b) we plot a(s) and a(s)r)(s) of which the latter is proportional to the fractional 
energy loss rate of the "leading particle", for ICS and TPP. Fig. 5 shows all the rates at 
redshift z = that affect the photons and electrons in the energy range we consider. 



3.1.5 Other Processes 

Other interactions that are neglected in this paper are all processes involving the productions 
of one or more e~e + pairs substituted by muon, tau lepton, and pion pairs, double Compton 
scattering (e7& — > 677), 77 scattering (77^ — > 77), Bethe-Heitler pair production (jX — ► 
Xe + e~ , where X can be an atom, an ion, or a free electron), the process 77^ — > e + e~7, and 
pair production on a magnetic field (7I? — > e~e + ). The total cross section of single muon pair 
production (77^ — > yU _ /i + ), for example, is smaller than electron pair production by about 
a factor of 10. Energy loss rates for TPP involving heavier pairs are suppressed by a factor 
~ (m/m e ) 1//2 in the limit of large s. Similarly, double pair production involving heavier pairs 
is also negligible |73J]. The double Compton scattering cross section is of order a 3 and must 
be treated together with the radiative corrections to ordinary Compton scattering of the 
same order. The corrections to the lowest order ICS cross section by processes involving m 7 
additional photons in the final state, e^b — > e + (m 7 + 1)7, m 7 > 1, turn out to be less than 
10% in the energy range under consideration f8~0fl . A similar remark applies to corrections 
to the lowest order PP cross section by the processes 77^ — > e + e~ + m 7 7, m 7 > 1. Photon- 
photon scattering can only play a role for z > 100 and energies below the redshift dependent 
pair production threshold Eq. ([]]) jSl], 82| . A similar remark applies to Bethe-Heitler pair 



production on atoms, ions and free electrons ||82|| . Pair production on a magnetic field of 
order 10~ 6 G which is typical for the field of Our Galaxy, is only relevant for E > 10 24 eV. 
This critical energy is even higher for the EGMF and this process is thus negligible in the 
analysis. 



3.2 Nucleons 

There are three major processes that affect the propagation of protons and neutrons: Electron- 
positron pair production by protons (PPP; p^t, — > pe~e + ), photopion production (A^, — ► 
N(mr),n > 1), and neutron /3-decay (n — > pe~v e ). 



3.2.1 Pair Production by Protons 

PPP provides the main energy attenuation for protons with energies below the GZK cut- 
off ||8"3"|. The energy threshold for this process is 

m e (m N + m e ) ^ u / j_n ~ 



E th = ^ ~ 4.8 x 10 i4 I — ) eV. (28) 
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Thus, for a microwave background photon (e ~ 10~ 3 eV), PPP ensues at a proton energy 
E ~ 5 x 10 17 eV. Below this energy, the protons cool essentially only by redshifting with the 
expansion. 

The PPP total cross section behaves very similarly to that for triplet pair production 



because PPP is almost identical to TPP and the expression for the total cross section 
away from the threshold may be given by Eq. (|24|). However, while TPP near its threshold 
is dominated by other processes, the exact behavior of PPP rates near the threshold are 
very important because PPP dominates the proton energy loss in that energy range. I use 
the parametric fits given in Ref. |85| for the cross section and inelasticity. Then I use the 



same approach in calculating the differential rates as we did for TPP. It can be shown that 
these rates are well approximated by a power law. On the other hand, the proton spectrum 
evolution due to PPP is well described by CEL because the inelasticities are smaller than 
10~ 3 at all relevant energies. Production of heavier pairs like n + fi~ is suppressed similarly 
to the case of TPP. The energy ranges for the produced pairs and the recoiling proton are 
given in Appendix A. 

3.2.2 Photopion Production 

Photopion production provides the main energy attenuation for nucleons above E ~ 10 19 eV. 
The energy threshold for this process is 

E th = mNm * + ml/2 * 6.8 x 10 16 (-V) ^ eV. (29) 

Thus, for a microwave background photon (e ~ 10~ 3 eV), photopion production ensues at a 
nucleon energy E ~ 7x 10 19 eV. Since publications on numerical studies of nucleon and 7-ray 
propagation usually do not contain detailed information on the implementation especially 
of multiple pion production, I present our approach here in some detail. First, I define a 
few suitable kinematic variables which depend only on the incoming particles. If E l ^ h is the 
photon energy in the laboratory frame (LF) where the nucleon is at rest, is the nucleon 
mass, and s is the squared center of mass (CM) energy, then the following relations hold: 

_ E^ b eE N 

K = — ^ = — r (l -fi) , 

m N my 

s = m 2 N {l + 2n). (30) 

Since laboratory measurements of cross sections are usually given in terms of E l ^ h , I will 
conveniently express everything in terms of k in the following. 

Concerning single pion production I consider the following reactions: 

1 + N - n° + N: ^ = E (*7 , (31) 

j=0 

1 + P - * + + n: ^ = (l~(3*xT 1/2 j:^)^y , (32) 
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7 + n ^ 7r-+p:^ = (l-/9V)- 1/a f;a v (ic)(i*) J '. (33) 

i=o 

The differential cross sections for these processes are expressed here in terms of k and the 
CM quantities Q*, (3*, and x* which denote solid angle, pion velocity, and the cosine of the 
scattering angle, respectively. The functions aij(n), i = 1,2,3 and j = 1, • • •, N are fitted 



to laboratory cross section data and I use fits up to order K = 3 ||86|| . The expressions in 
Eqs. (|3T|)-(|33|) can be easily rewritten in terms of the energies of the outgoing nucleons and 
pions in the cosmic ray frame (CRF) which I denote by E' a for a = p, n, tt + , tt~ , tt°. The 
relevant formulae are be given in Appendix B. 

Note that in Eq. ([31]) I have assumed identical cross sections for the two charge retention 
processes involving protons and neutrons. This is a very good approximation (see, e.g. 



Ref. |36|). Reactions ( p2[) and (^) constitute the charge exchange reactions for single pion 
production. 

We now turn to multiple pion production. Let us first consider the channel 7+p — > ir~+X 
where X stands for anything. This channel has been discussed in detail in Ref. |87fl . There, 
the Feynman x- variable x = Py/p^ax was introduced, which is the fraction of the pion parallel 
momentum p| in the CMF to its maximal value 

Pmax = m ^ ( 1+ ^)l/2 > ( 34 ) 

where D = [(k — e 2 /2) 2 — e 2 ] 1 ^ 2 , and e = m^jm^ is the ratio of the pion and nucleon 
mass. Denoting the transverse momentum with p± and the 7r~ energy in the CMF by E*- , 
the differential cross section for 7r~ production was written in terms of a structure function 

d 2 a lp ^ n - x = 7r-^ f(x,p 2 L ,s)dxdp 2 l . (35) 

Performing Lorentz transformations into the CRF where the proton and 7r~ energies are E p 
and E' n - (see Appendix B), this can be written as 



-> n -X 271171 N p 



7P— >7T — X 

d~EL. 




(36) 



1 + 2n) + 2p; ax j (1 + 2k)^x,s 

where x min > — 1 is chosen such that p\ > which is the second argument of /, and x x 
such that pi/feax) 2 + ^ 1- 



For s — >• oo the structure function / is independent of s |87|, |88fl . I will therefore neglect 
any s-dependence altogether. Furthermore, I take into account these processes only above 
some threshold which is sufficiently high such that the contribution of single pion production 
is negligible; I take k > 2(e + e 2 ). 
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Finally, for our purposes I assume that the remaining dependence of f(x,p\) factorizes 
into an x-dependent part and an exponential dependence on p\: 



f(x,p 2 ± ,s) 



A 2 



exp[~p 2 JA 2 ]f(x) 



(37) 



Here, A ~ GeV/6.4 is roughly of the order of the QCD scale and f(x) can be fitted to the 
data presented in Ref. |S7[ . 

Within these approximations we have finally 



dEL- 



da. 



■yp—*n X 




2~ f in x \ 2 ^ r-'''" 



m N 
A 



-e 2 + 



'EL 



E„ 



dxf(x) 



'1 + 2k) + 2D 



(38) 



'EL 



E n 



x 



where x min and x max are given by 



max 



x n 



_ s 2 -(E^/E p ) 2 (1 + 2k) 
2D(E^/E P ) 

,1 + 2k\ 1 ' 2 



E'„ 1 + 2k 
E p D 



At this point it is important to realize that f(x) = f c (x) + f p (x) 
a contribution f c (x) from the "central" pions and a contribution f p (x 



(39) 

can be divided into 
from production of 

multiple p mesons (sometimes also called the leading pion contribution) which subsequently 
decay into equal distributions of tt + and n~ . Therefore, f p (x) exclusively contributes to the 
production of 7r + and tt~ and corresponds to a charge retention process where the nature of 
the nucleon is unchanged. In contrast, f c {x) describes a process resulting in approximately 
equal distributions of tt° and tt~ with the probability for change of nucleon isospin being 
about 2/5 (from simple quark counting). 

From these assumptions it follows immediately that da^,/dE' n0 = da 7P ^ox/ 'dE' n0 is ob- 
tained by substituting f c (x) for f(x) in Eq. (08|) . 

In addition, I assume that inclusive and leading pion production takes place with the 
approximately constant cross sections <7 to t — 127 /A> and a p ~ 21.5 fib, respectively. We can 
then define the average central tt~ multiplicity by 



«->(«) 



7T 



O"tot 



dp]_ / dxf c (x 



exp[— p\/ A 2 
A 2 




o 1 + 2k 

x 2 + 



-1/2 



(40) 



The integration range is determined by pi/(p max ) 2 + x 2 < 1. By evaluating this formula 
one can see that the multiplicity increases asymptotically logarithmically in s for 

s — > oo. Applying charge conservation to the central pion distribution, making the above 
assumptions and in addition assuming tc + and tc~ distributions to be proportional to each 
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other uniquely determines da^/dE' iv+ = do lp ^ n +x/dE' n+ . It is obtained by substituting 
[1 + (2/5) (1 - o p /a tot )/ «_) («)] / c (x) + for /(x) in Eq. ©. 

It is now easy to compute the fractions r c (/c) and r p (/c) of the incoming nucleon energy 
which go into the central and leading pions, respectively (see Appendix B). Fig. 6 shows 
these fractions and the central and total 7r~ multiplicities as functions of s. Assuming a flat 
distribution for the outgoing nucleons, we then have 



da, 



dE' p 



da, 



p^n 



dEL 



On 



2r p {K)E p 







K 



+ - 



3 Otot - CTp q 



5 2r r (K)E„ 



2 a tot 



o n 



5 2rAK)E. 



6 



^-l+2r c ( 



^ - 1 + 2r c ( K ) 



(41) 



for the charge retention and charge exchange cross sections, respectively. For the processes 
involving an incoming neutron, I assume that the cross sections are also given by the above 
expressions after substituting p n and tt + <-> n~ everywhere. Fig. 7 shows the differential 
cross sections for production of 7r~, 7r + , 7r°, protons, and neutrons for an incoming proton 
for two different CM energies resulting from the formalism adopted above. Fig. 8 shows the 
inclusive pion production cross sections for nucleons as a function of s. Fig. 9 shows all the 
rates at redshift z = that affect the nucleons in the energy range we consider. 

Pions produced by nucleons quickly decay to EM particles and neutrinos and feed the EM 
cascade. 7r° decays into photons (tt° — * 77), and 7r ± decays to produce electrons, positrons, 



and neutrinos (n ± — > / u ± ^(z/ M ); y u ± — > e ± u e {v e )V ll {v^)) |89[]. Since the decay time of pions 
is very short compared to the timescale in the problem, I assume that pions are converted 
into secondary particles instantaneously. The decay spectra of the secondary particles may 
be calculated easily [pO[l. The expressions for the decay spectra are given in Appendix C. 



3.2.3 Neutron /3-decay 

Below ~ 10 20 eV, neutron /5-decay is the fastest process among the interactions that affect 
nucleons in the problem. The neutron decay rate is Y = r /7n = l/r n 7 n , where r n is the 
neutron lifetime (r n ~ 888.6 ± 3.5 sec), and 7„ is the neutron Lorentz factor. The range R n 
of a neutron is given as 

R n ~ I = cr nln ~ 0.9 (— !^ Mpc. (42) 



r V10 20 eV 

In calculating the spectrum of secondary particles, I neglect the proton kinetic energy in 
the neutron rest frame, as is usually done. The result can be found in standard textbooks 
such as Ref . [[Xj . 

Fig. 10 shows the energy attenuation lengths for cascade photons and nucleons as func- 
tions of energy in the CEL approximation. 



17 



4 Comparison with other Work 



Before I apply the propagation code to specific HECR injection scenarios in the next section, 
I compare the predicted spectra with results from other investigations for some standard 
situations. For a discrete source producing a differential injection spectrum F a (E) of particle 
type a (in units of number per energy per time) at redshift z = Z{, we obtain the spectrum 
j a {E) (in units of number per area per time per solid angle per energy) observed at z = in 
the following way: I impose the boundary condition 



N a (E, Zi 



-F a {E) 



(43) 



where r, is the comoving dimensionful source distance corresponding to redshift z% (r* = 
2Hq 1 \\ — (1 + Zi)" 1 / 2 ] in our cosmology), and solve the propagation equations for vanishing 
source terms. If we denote the resulting distribution at z = by N a (E), then j a (E) = 
N a (E) / (4k) and the modification factor M a (E, Zi), defined as in Ref. |42|, |4j], is given by 



M a (E, Zi 



F a (E) 



[l + Zi 



■ N a (E) 
N a (E, Zi ) 



(44) 



where I used the luminosity distance d L = r^l + z^ pTf . 

In Fig. 11 I plot the modification factors as defined in Ref. for discrete sources in- 
jecting protons with a power law at a given distance along with the corresponding curves 
51]. It can be seen that our results lie somewhat between results from Refs 



from Ref. 



51 



and j|7|] . In Fig. 12 I compare the nucleon, 7-ray and neutrino fluxes computed for monoen- 
ergetic proton injection at a given distance with results from Ref. |^TJ. In our prediction the 
secondary 7-ray flux at the low energy side is higher than the one given in Ref. |HJ by a 
factor ~ 10. I attribute that to the fact that the differential multiple pion production cross 
section used in our analysis (see Fig. 7) peaks at low energies. In Fig. 13 I consider the case 
of power law injection by a single source and compare the nucleon and 7-ray fluxes with cor- 
responding results in Ref. . The nucleon fluxes agree well, whereas, again, our prediction 
for the 7-ray flux is higher at the low energy side and lower at the high energy side. Since 
Refs. j|7| 51j do not give detailed information on their treatment of pion production, it is 
hard to give an exhaustive explanation of these differences. This, however, will not have an 
influence on our considerations where secondary 7-ray production by nucleons plays a minor 
role. 



5 Application to Models of HECR Origin 

We are now in a position to compute the cosmic and 7-ray fluxes predicted by various 
models of HECR origin. Since there is currently no unambiguous information on HECR 
composition, I will normalize the predicted sum of 7-ray and nucleon fluxes to the observed 
HECR flux. This is done to optimally enable an explanation for the events above 10 20 eV 
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without overshooting the UHE flux at lower energies (which might be explained by more 
conventional components) or predicting an excessive integral flux above 10 20 eV. I estimate 
the uncertainty in the predicted 7-ray flux at lower energies induced by this normalization 
procedure to be less than a factor ~ 3. 



5.1 GUT Scale Physics Models 

As already mentioned in the introduction, it has been suggested that HECRs may have 
a nonacceleration origin []27|, ^8], |2^, |3(J, [31], [32], |33|, [34| such as the decay of supermassive 



elementary "X" particles associated with Grand Unified Theories (GUTs), for example. 
These particles could be radiated from topological defects (TDs) formed in the early universe 
during phase transitions caused by spontaneous breaking of symmetries implemented in 
these GUTs (for a review on TDs, see [0]). This is because TDs, such as ordinary and 
superconducting cosmic strings, domain walls and magnetic monopoles, are topologically 
stable but nevertheless can release part of their energy in the form of these X particles due 
to physical processes like collapse or annihilation. The corresponding injection rate of X 
particles dnx /dt as a function of cosmic time t is usually parametrized as 

^ocT 4 ^, (45) 
dt 

where p > depends on the evolution of TDs. For example, X particle release from a 
network of ordinary cosmic strings in the scaling regime would correspond to p — 1 if one 
assumes that a constant fraction of the total energy in closed loops goes into X particles 



31]. Annihilation of magnetic monopoles and antimonopoles [27, p3| predicts p = 1 in the 



matter dominated and p = 3/2 in the radiation dominated era whereas the simplest 



models for superconducting cosmic strings lead to p = pqj . A constant comoving injection 
rate corresponds to p = 2 and p = 5/2 during the matter and radiation dominated era, 
respectively. 

The X particles with typical GUT scale masses mx of the order of 10 16 GeV subsequently 
decay into leptons and quarks. The strongly interacting quarks fragment into a jet of hadrons 
which results in mesons and baryons that are typically of the order of 10 4 — 10 5 . It is assumed 
that these hadrons then give rise to a substantial fraction of the HECR flux as well as a 
considerable neutrino flux. 

The shapes of the nucleon and 7-ray spectra predicted within such TD models are thus 
expected to be universal (i.e., independent of the specific process involving any specific kind 
of TD) at ultrahigh energies and to be dependent only on the physics of X particle decay. 
This is because at HECR energies nucleons and 7-rays have attenuation lengths in the cosmic 
microwave background (CMB) which are small compared to the Hubble scale. Cosmological 
evolutionary effects which depend on the specific TD model and are usually parametrized by 
Eq. ( f45|) are therefore negligible. In contrast, the predicted neutrino flux and the 7-ray flux 
below the pair production threshold on the CMB [see Eq. (|I])] depend on the energy release 
integrated over redshift and thus on the specific TD model. 
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I now discuss the particular form of the particle injection spectra expected from X particle 
release. I assume that each X-particle decays into a lepton and a quark each of an energy 
approximately half of the X particle mass mx- For reasonable extragalactic field strengths, 
the lepton (which I assume to be an electron in the following) will quickly be degraded 
by synchrotron loss producing synchrotron photons of a typical energy given by Eq. (|8|). 
This energy is typically much smaller than 10 20 eV where the resulting contribution to the 
7-ray flux is likely to be buried below the charged CR flux. For that reason, the GUT-scale 
lepton was usually omitted. However, for high EGMF strengths the synchrotron peak can 
approach 10 20 eV and thus could become relevant. For the present analysis I will thus include 
the source term for the GUT-scale lepton by writing its injection flux at energy E and time 
t as 

9 e (E,t) = ^^S(E-m x /2), (46) 

in units of particles per volume per time per energy. 

The quark from X particle decay hadronizes by jet fragmentation and produces nucleons, 
7-rays and neutrinos, the latter two from the decay of neutral and charged pions in the 
hadronic jets. The hadronic route is expected to produce the largest number of particles. 
The resulting effective injection spectrum for particle species a from the hadronic channel 
can be written as 

9ml E,t) = (47) 
at mx dx 

where x = 2E/mx, and dN a /dx is the effective fragmentation function describing the pro- 
duction of the particles of species a from the original quark. 

The spectra of the hadrons in a jet produced by the quark are, in principle, given by 
quantum chromodynamics (QCD). Suitably parametrized QCD motivated hadronic spectra 
that fit well the data in collider experiments in the GeV-TeV energies have been suggested 



in the literature [27J. The total hadronic fragmentation spectrum dN^/dx is taken to be of 
the form [27 

dN h (x) = f i|x- L5 (l - xf if x < x < 1 
dx 1 otherwise 

where the lower cutoff Xq is typically taken to correspond to a cut-off energy ~ 1 GeV. The 
spectrum Eq. (|8|) obeys energy conservation, /* dxx(dNh(x) / dx) = 1. Assuming a nucleon 
content of ~ 3% and the rest equally distributed among the three types of pions, we can 
write the fragmentation spectra as [p2l BO] 



^1 = (0.03)^, (49) 
dx dx 

dN n + dN n ~ dN n o /0.97\ dN h (x) 



dx dx dx V 3 / dx 

From the pion injection spectra one gets the resulting contribution to the injection spectra 
for 7-rays, electrons and neutrinos by applying the formulae in Appendix C. 
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Independently of the spectral shapes of the predicted nucleon and 7-ray fluxes, the ques- 
tion for the absolute normalization of the injection rates drix/dt in Eqs. (|45|), ([16|) and fl47|) 
arises. It has been shown, for example for cosmic strings |M ^TJ and annihilation of magnetic 
monopoles and antimonopoles f33|, that at least some TD models are capable of producing 
an observable HECR flux if reasonable parameters are adopted. For the purposes of this 
paper I will therefore not consider this issue and simply adopt the normalization procedure 
mentioned above. 

TD models of HECR origin are subject to a variety of constraints mostly of cosmological 
nature. These are mainly due to the comparatively substantial predicted energy injection 
at high redshift [see Eq. (|45|)]. Note that more conventional CR sources like galaxies start 
to inject energy only at a redshift of a few. Using an analytical approximation for the cas- 
cade spectrum below the pair production threshold on the CMB resulting from X particle 
injection, one can derive constraints from cascading nucleosynthesis and light element abun- 
dances, CMB distortions, and the measured 7-ray background fiE , 
region 



in the 100 MeV 



56] , as well as from observational limits on the 7-ray to charged CR flux ratio between 



10 13 eV and 10 14 eV |J3 
in Refs. 



5j|. The 100 MeV 7-ray background constraint was first discussed 



9], [55]. In the context of top-down models it was applied in Refs. [0, |95| on the 
basis of analytical approximations. 

In addition, there has been a claim recently |nj that TD models might be ruled out 
altogether due to overproduction of 7-rays in the range between the knee and ~ 10 19 eV. 
This would occur for EGMFs stronger than about 10~ 10 G due to synchrotron radiation 
from the electronic component of the TD induced flux which was normalized to the observed 
flux at 3 x 10 20 eV. However, in my opinion, the argument in Ref. |nj suffers from several 
shortcomings: First, monoenergetic injection of protons and 7-rays was used instead of the 
more realistic injection spectra such as the ones discussed above in Eqs. ( PD and (H). And 
second, only the case of a single, discrete TD source at a fixed distance from the observer 
was considered instead of more realistic source distributions and evolution histories. Finally, 
electron deflection due to the EGMF which can influence the processed spectrum from a 
single source was neglected. Nevertheless, I simulated the situation of Fig. 13 in Ref. ]51 
for an EGMF of 10 -9 G on which their claim is based. As a result I got a spectrum whose 
shape is roughly similar to Fig. 13 in Ref. [51], but the details of the spectrum differed 



somewhat, part of which can be attributed to a different model of the radiation background. 
The synchrotron peak I got was about an order of magnitude lower than in Ref. [51] relative 
to other parts of the spectrum. Most importantly, however, we observe that if the spectrum 
is normalized to the highest energy event this model would predict simply too many events 
above 10 20 eV including the original injection peak. Therefore, the model adopted in Ref. ]5l 



is not a realistic model for UHE CRs to start with. I thus conclude that it is not possible to 
rule out TD models on the basis of the discussion in Ref. 



51 



Our goal here is to reexamine the constraints based on the predicted 7-ray flux in the 
regimes around 100 MeV, between 10 13 eV and 10 14 eV, and between the knee and 10 19 eV, 
using our numerical techniques discussed in the previous section, I base this on realistic 
injection spectra and histories as discussed above. To my knowledge, this has not been done 
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yet despite its importance for making 7-ray flux based constraints more reliable. 

The redshift range of energy injection contributing to the 7-ray flux at energy E 1 today 
is given by 1 + z < (E th (z = 0)/E y ) 1/2 where E th (z = 0) is the PP threshold on the CMB 
at z = [see Eq. ([]])]. Since our interest is in the 7-ray flux at E 1 > 100 MeV, I maximally 
integrate up to 1 + z max = 10 3 . The spectrum in this energy range converges before we 
reach this redshift. A word of caution is in order for the predicted neutrino spectra. The 
UHE neutrinos interact with the universal neutrino background with T v ~ 1.95 K, and 
produce 11 where I = e, //, r, u, q, . . . via Z resonance |32], f4(J. The decay products of /1, r, q 
contain secondary neutrinos. Here I consider only simple absorption of UHE neutrinos, i.e. I 



integrate up to the average absorption redshift z^ due to this interaction |32|1 . The neutrino 
spectra also converge rather fast with increasing redshift for the parameters I used for TD 
models. Furthermore, the modification to the neutrino spectra due to the cascading by the 
aforementioned interaction is expected to be small for these parameters HO]. Therefore, the 



neutrino spectra given in this paper are expected to be good approximations to the real 
converged spectra. I leave a more detailed calculation of the UHE neutrino flux to future 
work. 

I performed simulations assuming uniform injection rates given by Eqs. for 
rrix = 10 23 eV and an injection history given by Eq. (|45|) for p = 1 (representative of scenarios 
based on ordinary cosmic strings and monopole-antimonopole annihilation) and for constant 
comoving injection (p = 2). Fig. 14 shows the results for a negligible EGMF and assuming 
our IR/O background model. Note that for a vanishing EGMF the 7-ray flux dominates the 
nucleon flux at UHEs and is higher by about an order of magnitude compared to predictions 
within the CEL approximation [compare Figs. 14(a) and 15]. This is due to the influence of 
non-leading particles on the development of the EM cascade. Fig. 16 shows the dependence 
of the results on the EGMF and the IR/O background. For an EGMF strength > 10~ n G, 
the 7-ray flux is determined by photon absorption and is thus harder. It is suppressed below 
a few 10 20 eV and dominates at higher energies which is in contrast to the case of a negligible 
EGMF [compare Figs. 14(a) and 16(a), (b)]. This scenario has the potential of explaining a 
possible gap in the HECR spectrum ||13|| . On the other hand, the neutrino flux is typically at 
least one order of magnitude larger than the other components. However, we note that the 
probability that these UHE neutrinos generate a shower in the atmosphere is smaller than 
10~ 5 0]. We also find that the predicted integral neutrino flux above ~ 10 20 eV is about 
10 4 times smaller than current limits from the Frejus experiment |78| . 

In the case of a negligible EGMF and absence of an IR/O background [see Fig. 16(c)], our 
normalization procedure leads to a 7-ray background below ~ 10 14 eV which is about 20 times 
lower than analytical estimates adopting a normalization based on the CEL approximation 
for the 7-ray component alone |5B[]. This is caused by the aforementioned influence of 
the non-leading EM particles on the UHE flux on which the normalization depends. For 
given backgrounds, EGMF strength, and flux normalization at HECR energies, the 7-ray 
background flux below ~ 10 14 eV is proportional to the total energy injection which increases 
monotonically with decreasing p. Comparison with the 7-ray background observed around 
100 MeV |5D| and recently up to ~ lOGeV jES), |6(J clearly rules out the cases with p < 1 
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within our IR/O background model and negligible EGMF [see Fig. 14(a)]. For an EGMF 
near its currently believed upper limit ~ 10~ 9 G | 79fl , proper normalization of the different 
predicted spectral shape at UHEs leads to an increase of the predicted low energy 7-ray 
background by about a factor 5, thus tightening the constraint somewhat (Fig. 16). The 
7-ray flux level between ~ 10 11 eV and ~ 10 14 eV is very sensitive to the IR/O background, 
and in the extreme case of absence of any IR/O flux it increases by about a factor 10 2 relative 
to the level predicted by our IR/O background model. At the same time, the flux below 
~ 10 GeV goes down by about a factor of 10 for vanishing IR/O flux [see Fig. 16(c)]. I stress 
that in any case, the scenarios considered here are currently neither constrained by the limit 
on the 7 to charged CR flux ratio below 100 TeV |[93|| , nor by the synchrotron peak between 
the knee and ~ 10 20 eV. 

Analytical arguments ]56| suggests that for a given normalization of the spectra at the 
highest energies and a given injection history, the total injected energy and thus the 7-ray 
flux below the PP threshold is roughly proportional to m x ~ i '. Here, a fragmentation function 
N a (x) was assumed which is roughly proportional to x~ q for x > E ohs /m x , i.e. q = 1.5 in 
the case of Eq. (^). This allows one to rescale the above constraints to different values of 
mx- For example, for p = 1 the constraint on q is roughly 

q ~ 2 ~ 3 + log (ml/10 23 eV) ' (50) 

However, the effect of the cascading and the EGMF complicate the problem considerably 
because the UHE spectrum depends sensitively on those effects. More accurate estimates 
can be achieved only by a separate numerical simulation of the case mx 3> 10 23 eV. I leave 
that to a forthcoming letter which will summarize the results [p6 . 

In order to mitigate or avoid overproduction of the 7-ray background, based on analytical 
considerations of the cascade spectrum, Chi et al. |95| have recently suggested somewhat 
different injection spectra. In this case injection of nucleons, 7-rays, and neutrinos is again 
given by Eq. ([|7]) where the following fragmentation functions are adopted: 

dN N (x) _ L5 dN v {x) _ dN^x) _ 2A 

; — Aw x j ; — ; — ^7 x > 

dx dx dx 

with 

±* 0.028 (iO^eVr. (52) 
A N \ m x J 

The condition (|52] ) comes from the requirement [£5| that the photon-to-nucleon (7/iV) ratio 
at injection at energy E = 10 20 eV, i.e., at x = 2 x 10 20 eV/mx be ~ 60. The spectra 
Eq. ( |STD are absolutely normalized such that the total quark energy mx/2 is injected between 
E = 5 x 10 19 eV and and E = m x /2. 

Fig. 17 shows results obtained by assuming the fragmentation functions given by Eqs. (|5T|), 
(|52]). The 100 MeV —10 GeV 7-ray background constraint is basically unchanged from the 
case of the QCD motivated injection spectra for vanishing EGMF on which it depends more 
weakly. Note that this scenario has the potential to explain a HECR spectrum continuing 
beyond 10 20 eV without any break or gap |I3]| . 
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5.2 Gamma Ray Burst Models 



Recently, it has been suggested that UHE CR could be associated with cosmological GRBs |24], 
|2~5| , |26|| . This was mainly motivated by an apparent numerical coincidence: Assuming that 
each (cosmological) GRB releases an amount of energy in the form of UHE CRs which is 
comparable to the total 7-ray output normalized to the observed GRB rate (about 10 51 erg 
per burst), the predicted and the observed UHE CR flux at the Earth are comparable. It 
should be mentioned, however, that it is not clear whether constraints on cosmological GRB 



distributions are consistent with HECR observations [|9j] . 

In these models protons are accelerated to UHE via first order Fermi acceleration. Since 
there are no firm predictions for the injection spectrum, I assume the hardest possible spec- 
trum proportional to E~ 2 up to a maximal energy of 10 23 eV. Furthermore, assuming a 
constant comoving injection rate up to some maximal redshift z max , we can write 

%(E, t) oc r 2 E- 2 Q(z max - ^)6(10 23 eV - E) . (53) 



The authors of Ref. [p6j pointed out that bursting sources in combination with deflection of 
protons in the EGMF could lead to UHE CR spectra with a time variability on a scale of 
~ 50 yr. This might allow reasonable fits to the observed HECR spectrum. However, I only 
consider the continuous injection of CRs in this paper for illustrative purposes. Since the 7- 
ray background depends only on the average flux, the only uncertainty in its flux level comes 
from the fit to the HECR events. I estimate the uncertainty introduced by normalizing the 
average flux to be less than a factor 1.5. Fig. 18 shows the results for various values of z max . 
I conclude that these models are currently unconstrained by the 7-ray background, although 
they still have difficulty explaining the highest energy CR events. 



6 Conclusions 

I have performed detailed numerical simulations for the propagation of extragalactic nucle- 
ons, 7-rays, and electrons in the energy range between 10 8 eV and 10 23 eV. My goal was to 
explore constraints on various models of HECR origin from a comparison of predicted and 
observed 7-ray fluxes at lower energies. The main focus thereby is on models which associate 
HECRs with GUT scale physics or with cosmological GRBs. 

I find that at present the TD scenarios are primarily constrained by the observed 7-ray 
background between ~ 100 MeV and ~ 10 GeV but not by the limit on the 7 to charged 
CR flux ratio below 100 TeV. The CEL approximation usually does not take the IR/O 
background into account, and thus may not be directly compared to the numerical calculation 
because the presence of the IR/O background may affect the 7-ray flux level at 100 MeV 
by an order of magnitude. There is also a significant difference for the UHE spectrum 
between predictions by the CEL approximation and my numerical simulation. For an EGMF 
strength < 10~ n G the TD models yield the 7-ray flux which is at about the same level as 
or below the current observed flux, depending on the adopted parameters. On the other 
hand, an EGMF stronger than ~ 10~ n G stops the cascade at UHEs and the UHE end of 
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the spectrum is suppressed significantly. Thus, the level of the 7-ray flux at about 100 MeV 
is higher relatively, tightening the constraints. However, these results are rather insensitive 
to different models of the IR/O background [63], although they are somewhat dependent on 
the poorly known universal radio background flux. 

I conclude that TD scenarios with QCD motivated injection spectra up to energies < 
10 23 eV are still viable if injection occurs uniformly or from a discrete source. This is in 
contrast to a recent claim in the literature [|5l|. In case of uniform injection this assumes an 
injection history motivated by energy release from a network of cosmic strings in the scaling 
regime or from monopole-antimonopole annihilation (p = 1). Higher injection energy cutoffs 
are allowed for either a weaker source evolution or for injection spectra somewhat steeper 
than the QCD motivated spectra. For EGMF strengths larger than ~ 10~ 10 G, some of the 
predicted TD spectra have the potential to explain a possible gap in the HECR spectrum. 
The cosmological GRB scenarios recently suggested in the literature [24, 25. 26] are currently 
unconstrained by these limits. 

With the arrival of the anticipated Pierre Auger Cosmic Ray Observatories [57], it is 
expected that the UHE end of the CR spectrum will be known with much better accuracy. 
Constraints derived from the influence of CR propagation on the observed spectrum will 
then be one of the most powerful tools in discriminating between models of HECR origin. 
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Appendix A: Triplet pair production 



The expressions for the differential spectra of produced pairs and the recoiling electron 
(positron) for TPP by a very energetic electron on a soft photon can be found in many 
papers J75], [7|, [FT]]. I adopt the analytic approach used in [TfJ. For an interaction of an 



electron of energy E with a photon of energy e (E ^> e), the double differential cross section 
with respect to the positron produced with energy E + at solid angle fl + can be expressed as 

" a T ■ TTTT* , I / Ac^V, (54) 



dE + dVt p+ 167T 3 p ■ p( 

where p • k is the scalar product of the initial electron and photon four-momenta, p + is the 
magnitude of the produced positron three-momentum, Q p * is the solid angle of the recoiling 
electron, and pt and A t are given in Ref. ||76|| . The single differential cross section with respect 



to the positron energy may be obtained by integrating Eq. (|54|) over the positron solid angle 
Q + numerically. The differential cross section for the produced electron is identical to that 
for the positron due to symmetry. In doing the integral, it is useful to use the approximation 
where the dependence of the cross section on the azimuthal angle of the outgoing particles 
in the cosmic ray frame is neglected, as was mentioned before. 

Finally, we can obtain the total TPP cross section by integrating Eq. ( |54"|) numerically 
over the kinematic range of electron and positron energies given by 

p _ E tot (s ~ ml) ± P tot [s (s - Aml)}^ 

2s + mj 

where So = e(E + p), and E tot and P tot are the total incident energy and momentum, 
respectively. 

I also give here the kinematic energy range for the outgoing electron and positron and 
the recoiling proton in case of pair production by protons (Section 3.2.1): 

p _ gtot(go ~ m e X) ± P to t[gQ - 2m e g (m e + A) + m 2 e (m 2 e - m 2 N )/2} 1/2 

2s + m 2 N 

where so is as defined previously, A = [(m 2 , + m 2 N )/2] 1 / 2 . 



Appendix B: Photopion production 

First, I express the differential cross sections for single pion production in terms of the CRF 
energies E' a , where a = N,n: 



doi doi 2tt 1 + 2k 



k 



dE' N ' dE' n E N D — 



x 



1 1 for % — 1 (charge retention) 

(1 — f3*x*) for i = 2,3 (charge exchange) 
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where /3* = D/(k + e 2 /2) and x* can be expressed as a function of E' N or E' n and i? 



AT- 



X 



'1 + 2k){E' n /E n )-\- K-e 2 /2 
D 



[l + 2K)(E'jE N )-K-e 2 /2 
D 



(58) 



Finally, I compute the fractions t c (k) and t p {k) of the incoming nucleon energy going into 
the central and leading pions. These fractions are given by integrating the differential cross 
section da/dE n for the respective process, weighted by the pion energy E n , in the CRF, and 
dividing by the corresponding total cross section. Using Eq. fl35|), 



and Eq. 
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), we end up with 
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(59) 



(60) 



-1/2' 



Again, the integration ranges are obtained by the requirement p\j '{p* m&x ) 2 + x 2 < 1. The 
formula for r c (/c) can be obtained from this by substituting o p — > cr tot — o p and 2j p {x) — > 
[3 + (2/5)(l - (T p /a %ot )/ (n%-) (k)] / c (x), where («) was given in Eq. 



Appendix C: Pion decay spectra 

First, we define the decay spectrum N a (E) as the differential number of the secondary particle 
a at energy E. Then the spectrum is normalized as 

J dEN a (E) = n a , (61) 

where n a is the number of particles a produced by decay of a single pion. For example, 
n 7 = 2 for 7T° decay. 

First, the photon spectrum from 7r° decay is 

N y (E) = J- for E < E n . (62) 
E n 

Before we calculate the charged pion decay spectra, we note the fact that the pions pro- 
duced from photopion production are always relativistic. Thus, we may make the relativistic 
approximation for both pions and the resulting muons. For example, the decay spectra of 



32 



7r + read: 



N + (E) ~ N- (E) ~ { ^ {A ° + + A3Z3) f ° r E ~ 

eH } iV ^- j _J_ {BQ + B > o \ nz + B2Z 2 + B ^ ) for rE n <E< E„, 

N^(E) = -L^- hrE<(l-r)E w (63) 



( TZ ^ 17 (C + C 2 ^ 2 + C3Z 3 ) for E < rE w 

Tr ^ K (D + D' \nz + D lZ + L> 2 z 2 + D 3 z 3 ) for rE n < E < E w , 



where r = m 2 /m 2 , z = E/E n , and coefficients are given as 



(A , A 2 , A 3 ) = (0.94486, -2.7892, 1.2397), 

(B , B' , B 2 , B 3 ) = (-2.4126,-2.8951,4.3426,-1.9300), 
(C ,C 2 ,C 3 ) = (1.1053,-4.46883,3.71887), 

and 

(D , D' , £>!, D 2 , D 3 ) = (13.846, 5.37053, -28.1116, 20.0558, -5.7902). 

The average energies of the secondary particles are (E e +) = (Ep ) = 0.265^, (-E^) = 
0.257 E n , and (E U/x ) = 0.213E n respectively. The decay spectra from n~ are obtained by 
substituting particles accordingly. 
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e (eV) 



Figure 1: The universal background radiation intensity spectrum at z = (solid line) used 
in our model. The separate contributions from the radio (short dashed line), the IR/O (long 
dashed line) background, and the CMB (dotted line) are also shown. 
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Figure 2: The effective comoving density of radio and IR/O sources whose luminosities are 
normalized at z = 0, as a function of redshift. This corresponds to N c (z) in Eq. (|3|). The 
IR/O source density is assumed to cut off at z = 5. 
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Figure 3: Gyroradii (dashed lines) and synchrotron loss lengths (solid lines) of electrons for 
various strengths of the EGMF in units of gauss (G) as indicated. 
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Figure 4: The total cross sections, cr(s), and the cross sections times the average inelasiticity, 
a(s)r)(s), which is proportional to the fractional energy loss rate of the leading particle: (a) 
For PP (solid line and short dashed line, respectively) and DPP (dotted line). 
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Figure 4: (b) For ICS (solid line and short dashed line, respectively) and TPP (dotted line 
and long dashed line, respectively). 
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Figure 5: The relevant interaction rates at z ^rte^Wect the photons and electrons in the 
energy range we consider. The key is identical to the key for Fig. 4. The rates are calculated 
by folding the total cross sections and inelasticity weighted cross sections with the present 
background photon spectrum shown in Fig 2. 
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Figure 5: (b). 
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Figure 6: (a) The fractions r c and r p of the incoming nucleon energy which goes into the 
central and leading pions, respectively, as functions of s. 
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Figure 6: (b) The average central and total tt multiplicities, (n°_) [see Eq. (|40|)1 and (n n -) 
[resulting by substituting f c (x) — > /(x) in Eq. fl40D l as functions of s. 
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Figure 7: The differential cross sections for production of rc~ (dotted lines), n + (short dashed 
lines), 7T° (solid lines), protons (long dashed lines), and neutrons (dash-dotted lines) for the 
collision of a proton of energy E with a background photon at squared CM energy s, from 
the formalism adopted in Section 3.2.2: (a) For E = 3 x 10 20 eV, s = 2.1 GeV 2 . 
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Figure 8: The inclusive multiple pion production cross section for protons (solid line) and 
neutrons (dashed line) as a function of s. 
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Figure 9: The interaction rates and energy attenuation rates for multiple pion production 
(solid line and short dashed line, respectively) and PPP (dotted line and long dashed line, 
respectively). The rates were obtained by folding the cross sections and inelasticity weighted 
cross sections with the present background photon spectrum shown in Fig 2. 



46 




Figure 10: The energy attenuation lengths for cascade photons and for protons as a function 
of energy assuming the radiation background photon spectrum shown in Fig. 2. These 
curves were obtained by running the code over small distances and ignoring the production 
of non-leading particles, which corresponds to the CEL approximation. 
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Figure 11: The modification factors as defined in Ref. |44| for discrete sources injecting a 
E~ 2 proton spectrum extending up to 3 x 10 20 eV at a given distance d or redshift z resulting 
from our analysis (solid lines). Also shown are the corresponding curves from Protheroe and 
Johnson |5lJ (dashed lines): (a) For d = 256 Mpc. For further comparison with results from 
Refs. f|7|] see the discussion in Ref. ]5T]] . 
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Figure 11: (b) For z = 0.6. 
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Figure 11: (c) For z — 1. 
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Figure 12: The differential ffuxes of 7-rays (solid line), nucleons (long dashed line), neutrons 
(short dashed line) and u^, u^, u e , v e (thin solid lines in decreasing order) for monoenergetic 
proton injection at an energy E = 10 21,5 eV and a distance d = 32 Mpc: (a) Result from our 
analysis in arbitrary units; (b) Corresponding results from Ref. [^l[ . 
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Figure 13: The differential ffuxes of 7-rays (solid lines) and nucleons (long dashed lines) from 
a discrete source injecting a E~ 2 proton spectrum extending up to 10 22 eV, located at the 
redshift indicated: (a) Result from our analysis in arbitrary units; (b) Corresponding results 
from Ref. H7 
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Figure 14: Predictions for the differejafial flux>es'bf j-yuys (solid line), prot6ns 
line), neutr4© fs^rt dashed line)^and Vp(v\v e (yQ (thin solid line\in decreski 



injection rates ^^^fflffl,^,^ MMMM 

functions Ec, 

(a) p = 1. Also shown are the combined data from the My s Eye M,y\\ and the ALrAb/ 



g dashed 
g order) 



by a typical topological defect scenario for a vanishing EGMF. This model assumes uniform 

entation 
2® for: 

x i) p = 1. Also shWn^re^e'c^mbmecfclata from the Fly's Eye |^[| "and "the AG ASA. |TT|, |T2 
experiments above 10 19 eV (dots with error blfrs^fprecewise power law fits to the charged 
CR flux (thick solid line) and observational upper limits on the 7-ray flux around 100 MeV 
from Refs. |58|, |59|, |6(J (dotted lines in decreasing order). The arrows indicate the limits on 
the 7-ray to charged CR flux ratio from Ref. ||93||. 
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Figure 14: (b) p = 2. 
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Figure 15: Same as Fig. 14(a), but neglecting non-leading particles in the EM cascade 
development. Only the UHE part is important in comparison with Fig. 14(a). 
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Figure 16: (a) Same as Fig. 14(a), but for an EGMF of KT 9 G. 
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Figure 16: (b) Same as Fig. 14(a), but for an EGMF of 10" n G. 
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Figure 16: (c) Same as Fig. 14(a), but assuming absence of any IR/0 background. These 
panels demonstrate the EGMF and IR/O background dependence of the predicted HECR 
flux shape and composition and its influence on the prediction for the 7-ray background at 
lower energies. 
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Figure 17: (a) Same as Fig. 14(a), but assuming the fragmentation functions given by 
Eqs. (|51]), (11) with m x = 10 23 eV. 
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Figure 18: Predictions for the differential fluxefe of 7-rays (solid line), protons (long dashed 
line), neutrons (short dashed line), and u e , v e (thin solid lines from the top) by the 

GRB injection scenario given by Eq. ([53]) for vanishing EGMF for: (a) z max = 1. Observa- 
tional data and constraints are presented as in Fig. 14. 
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Figure 18: (b) z max = 4. 
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